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NONSYMMETRIC BALLASTIC RANGE, HEIGHT, TIME-OF -FLIGHT 
AND OPTIMAL FLIGHT PATH ANGLE COMPUTATIONS 
WITH PROGRAMS FOR A HEWLETT-PACKARD 65 CALCULATOR 



I . The Introduction and Apology . 

The range of a ballistic missile over the surface of the 
Earth is a topic which is discussed in almost every text on 
astronautics (Reference 1, for example). Unfortunately, none 
of the available texts contained a needed procedure for the 
exact solution of the unsymmetric ballistic problem (defined 
in Section II) ; all procedures found were approximate and 
were based on a symmetric approximation [See References 2 and 
3, for example] . The author of this report faced the dilemma 
of either (1) taking the time to solve the problem afresh, or 
(2) taking an inordinate longer time to do a thorough search 
to find a solution which might not be the one desired. The 
iormer course of action was taken and all due apologies are 
hereby extended to those that have published more elegant 
solutions . 

All of the computational methods are outlined in Section 
III. Instructions and a sample problem for the Hewlett- 
Packard 65 Calculator (HP-65) are given in Section IV. This 
HP-65 programs are in Section V, and the development of the 
3 opt procedure is in the Appendix. Familiarity with 
Newtonian two-body theory is assumed. 



II. The Problem 



Basically, the problem is to compute the range, S, and 
time-of-f light , t^, of a ballistic missile given the booster 
cutoff at height h„ after a vertical ascent, a cutoff 
velocity v L , and a subsequent flight path angle 8 (Fig. 1) . 
Given that a solution to the stated problem is available, the 
next natural question to ask is, "What flight path angle ® 0 p t 
will obtain the maximum range S ?" The solution is fairly 
simple in the symmetric ballistic problem, i.e. when h £ = 0, 
and called symmetric because the launch point L and target 
point T are symmetric with the line of apses of the 
Keplerian orbit (i.e. = (f >2 in Fig. 1). The unsymmetric 

ballistic problem is that in which h_ * 0 . Here we assume 

£j 

that h £ is given, the Earth with radius r E is spherical, 
non-rotating, and has no atmosphere. The effect of the Earth's 
rotation may be calculated figuring how far the target has 
moved during the time of flight of the missile. Explicitly, 
we assume that Keplerian orbit is obtained at height h„ with 

£j 

terminal booster cutoff velocity v_ and flight path angle 

1j 

8 (or 8^) , and the target T is at the Earth's surface. 

The computation of range and time-of-f light for given h E , 
r E' V L ' an< ^ ^ can easil Y b e obtained from the equations of 
elliptical orbits. The determination of 8 0 p t is somewhat 
complicated, however, and requires an iterative procedure. In 
all of the developments we assume that L and T are sepa- 
rated by the line of apses (which implies that 0 < 8 ^ • 
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Line of apses 




FIGURE 1 
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III. The Computational Procedures . 



A. Range 



1. Let r T and r m denote the distance 

Li 1 

of the launch point L and target point 
T from the center of the Earth. Usually 
we let r L = r E +h E and r T = r E 
although r L and r T need not be so 
defined in general. 

2. From r T and v T , compute the semimajor 

Li 

axis of the ellipse. 



gravitational constant for the Earth. We 
have assumed that the mass of the missile 
is negligible with respect to the mass of 
• the Earth. 

3. Using the flight path angle 3, compute 
the elliptical semiparameter 



a 




2 

where y = k £ and k £ is the (Gaussian) 




4. Compute the eccentricity 




4 



5. Compute the true anomaly of the launch 



point 




where 0 < f T < tt . Note that the line 
of apses can be said to have true anomaly 

f = IT . 

6 . Compute the true anomaly of the target point 




7. The surface range of the missile is then 

S = r E (f T ‘ f L> ' 



B. Time-of-Flight 

1. Compute the eccentric anomaly of the 
launch point. 




where 0 < E < tt . 

2 . Compute the eccentric anomaly of the 
target point, 




where tt < E T < 2ir . 
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3. Using Keplers' equation, compute the 



time-of -flight 



= 



a 3/2 

/y" 



(E t - E L )-e(sinE T - sinEjJ 



C. 



Height of Trajectory & Circular Velocity 

1. The maximum height of the trajectory above 
the Earth's surface for any h E , v L , and 
3 is given by 



>W = h = a(l + e) - r E 



2. Circular Velocity 

The velocity required for a circular 

orbit at r is 
L 



v o = Z 57 *] 



v will be used in Section III D, 
c 



D. Optimum Flight Path Angle & Maximum Range 

For a specified h and v T , a value of 

hi Li 

3 can be found such that 



S max - S<S opt> ’ S(6) ' 



In this unconstrained optimization problem 



we set 



as 

93 



3 = 3 



= 0 



and solve 



opt 
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for • The relationship between S and 

$ is given by the computational formulas for 
the range in part A. Although an explicit 
formula for B is unobtainable, the following 
exact relationship is useful: 




where a = — . 

c 



In the symmetric ballistic problem r L = r T ; 
the equation for B 0 p t then reduces to 



( 2 ) 



cos^B 



opt 



2-a 



2 ' 



which is the result given in References 1, 2, and 3. 



find 



In the unsymmetric case Equation (1) is used to 

B . as follows: 
opt 



1. Solve (2) for B . . 

opt 

2. Substitute B . into the right hand side 

opt 

of (1) to obtain a new value for B . • 

opt 
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3. Average the last two values of B Qpt to 

obtain a new value of and return 

opt 

this value of B 0 p t to step 2. 

4. It has been found that five iterations are 
sufficient to obtain a reliable value for 

^opt 

5. Using B compute the maximum range 

S max = S ^opt ) using procedure A. 
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IV. The Illustrative Example. 



The HP-65 program consists of four magnetic cards which 
are used as follows: 

Card 1: This card must be entered first. It is vised 

to set the two constants r_ and y for any 
choice of four distance units and any choice of 
four time units. The distance units are 
kilometers, meters, statute miles, and nautical 
miles; the distance unit is entered by pressing 
key A, B, C, or D, respectively. The display 
shows the equatorial Earth radius r £ in the 
respective unit. The time units are seconds, 
minutes, hours, and days; the time unit is 
entered by pressing key A, B, C, or D, respec- 
tively. The display shows the constant y in 
the appropriate (distance) /(time) unit. Note - 
after entering program Card 1, the distance unit 
must be entered before the time unit is entered 
since the keys A through D perform a dual func- 
tion. The basic definitions used are: 

1 statute mile = 1.609344 kilometers (km) 

1 nautical mile = 1.852 km 
r = 6378.160 km 

£j 

y = 398603 km 3 /sec 2 

Card 2: This card computes S, h , and v given 

Illd X O 

h^, v T and 8. It must be used after the 
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execution of Card 1 and it may be used before 
and/or after the execution of Cards 3 and 4. If 
it is used after Card 4 to compute s max ' then 
S 0 p£ as computed by Card 4 need not be re-entered. 
Otherwise, in any order: enter h„ and press A, 

£j 

enter v T and press B, enter 8 and press C. To 
find the range S, press E. Then to find ^ max ' 
press R/S; and then to find v , press R/S again. 
Note: If h„, v T , and 8 have been entered 
using Card 3 , then they do not have to be 
reentered when using Card 2. Also, any new value 
of h E , V L , or 8 can be entered and the range 
computed without having to re-enter unchanged 
values. 

Card 3: This card computes the time of flight, t^ . 

It must be used after the execution of Card 1 
and it may be used before and/or after the 
execution of Cards 2 and 4. If it is used 
after Card 4 , then ^ 0 p t as computed by 
Card 4 need not be re-entered; if it is used 
after Card 2, then only new quantities need be 
entered. Otherwise, in any order: enter h E 

and press A, enter v L and press B, enter 8 
and press C. To find the time of flight, press 
E. 
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Card 4: This card computes B q . It must be used 

after h^ and have been entered from 

either Card 2 or Card 3. No entries are 
required. To compute 8 press any key A 

through E; five iterations will be performed 
and the resulting value of ^ Q p t will be dis- 
played. To perform five more iterations, press 
R/S. To monitor the intermediate iterations 
press keys f SF2 . Each intermediate result will 
be displayed. To continue, press R/S. To dis- 
able the monitoring, press keys f ^ SF2 . 

Using the four program cards as described above find 
the range and time-of-f light for a ballistic missile that 
has a velocity of 10000 knots at a cut-off altitude of 50 
nautical miles for a flight path angle of 30° and 40°. 
Also, find the maximum possible range and time-of-f light 
for the same missile. 



Proceed as follows: 



Instruction 


Data 


Key 


Display 


Enter Card 1 








Distance unit is n.mi. 




D 


r E =3443.93n.mi 


Time unit is hours 




C 


y=8.13247xl0 11 


Enter Card 2 








Enter cutoff altitude 


50 


A 


r=3493.93n.mi 

U 


Enter cutoff velocity 


10000 


B 




Enter flight path angle 


30 


C 




Compute range 




E 


S=1928 . 08 n.mi . 
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Step 


Instruction 


Data 


Key 


8a. 


Optional: Display h m!iw 

max 




R/S 


8b 


Optional: Display v 

c 




R/S 


9 


Enter flight path angle 


40 


C 


10 


Compute range 




E 


11 


Enter Card 3 






12 


Compute time-of-f light 




E 




Note: hg = 50, v L = 10000, and 








6 = 40° since no further 








entries have been made. 






13 


Enter flight path angle 


30 


C 


14 


Compute time-of-flight 




E 


15 


Enter Card 4 






16 


Compute 6 opt (5 iterations) 




E 


17 


Make 5 more iterations to be sure 




R/S 


18 


Enter Card 2 






19 


Compute 

max 




E 


20 


Enter Card 3 






21 


Compute time-of-flight for 




E 




maximum range. 







Display 

h =344 . 74n .mi 
max 

v =15256. 47knot 
c 

S=1956 . 60n.mi. 
t^=0 . 3043hrs . 



0 . 2504hrs . 



8 opt =36 ° 0895 

S opt =36 ° 0895 

S ma X =1973 - 95n '" 



W 1 ' 0 ' 2846 



1 
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V. The Program Listings . 
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CARD 1 



CARD 2 



Set units of distance and time 
display r_ and p 

Ei 



Compute range, h^^, v c 



10 



20 



30 



40 



Key 


Code 


Key 


Code 


Key 


Code 


Key 


Code 


Entry 


Shown 




Entry 


Shown 




Entry 


Shown 




Entry 


Shown 


LBL 


23 




3 


03 




RCL 4 


3404 




X 


71 


B 


12 




7 


07 




+ 


61 




1 


01 


EEX 


43 




8 


08 




STO 3 


3303 




- 


51 


3 


03 




• 


83 




R/S 


84 




CHS 


42 


CHS 


42 




1 


01 




LBL 


23 




f 


31 


GTO 


22 




6 


06 




B 


12 




y 


09 


0 


00 




X 


71 




STO 2 


3302 




STO 6 


3306 


LBL 


23 




STO 4 


3304 




R/S 


84 




RCL 8 


3408 


C 


13 




R/S 


84 




LBL 


23 




RCL 3 


3403 


1 


01 


60 


LBL 


23 


10 


C 


13 


60 


D 


14 


* 


83 




B 


12 




STO 7 


3307 




RCL 8 


3408 


6 


06 




6 


06 




R/S 


84 




RCL 4 


3404 


0 


00 




0 


00 




LBL 


23 




D 


14 


9 


09 




GTO 


22 




D 


14 




+ 


61 


3 


03 




1 


01 




g 


35 




g 


35 


4 


04 




LBL 


23 




RAD 


42 




IT 


02 


4 


04 




C 


13 






81 




ENTER 


41 


GTO 


22 




3 


03 




1 


01 




+ 


61 


0 


00 




6 


06 




- 


51 




g x<— >y 


3507 


LBL 


23 


70 


0 


00 


20 


RCL 6 


3406 


70 




51 


D 


14 




0 


00 






81 




RCL 4 


3404 


1 


01 




GTO 


22 




f-1 


32 




X 


71 


• 


83 




1 


01 




COS 


05 




R/S 


84 


8 


08 




LBL 


23 




RTN 


24 




RCL 6 


3406 


5 


05 




D 


14 




LBL 


23 




1 


01 


2 


02 




' 8 


08 




E 


15 




+ 


61 


GTO 


22 




6 


06 




g 


35 




RCL 5 


3405 


0 


00 




4 


04 




DEG 


41 




f 


81 


LBL 


23 




0 


00 




RCL 7 


3407 




RCL 4 


3404 


A 


11 


80 


0 


00 


30 


f 


31 


80 


_ 


51 


1 


01 




GTO 


22 




COS 


05 




R/S 


84 


LBL 


23 




1 


01 




RCL 2 


3402 




RCL 1 


3401 


0 


00 




LBL 


23 




X 


71 




RCL 3 


3403 


g 


35 




A 


11 




RCL 3 


3403 




i 


81 


1/x 


04 




1 


01 




X 


71 




f 


31 


ENTER 


41 




LBL 


23 




ENTER 


41 




. 


09 


ENTER 


41 




1 


01 










* 




ENTER 


41 




ENTER 


41 




X 


71 




R/S 


84 


X 


71 




X 


71 




RCL 1 


3401 




GTO 


22 


X 


71 


90 


RCL 1 


3401 






81 




E 


15 


3 


03 




X 


71 


40 


STO 8 


3308 


90 


g NOP 


3501 


9 


09 




STO 1 


3301 




2 


02 




g NOP 


3501 


8 


08 




R/S 


84 




RCL 3 


3403 




11 


3501 


6 


06 




g NOP 


3501 




T 


81 




II 


3501 


0 


00 




g NOP 


3501 




RCL 2 


3402 




tl 


3501 


3 


03 




11 


3501 




ENTER 


41 




It 


3501 


X 


71 




11 


3501 




X 


71 




11 


3501 


STO 1 


3301 




11 


3501 




RCL 1 


3401 




11 


3501 


CL x 


44 




11 


3501 




V 


81 




11 


3501 


6 


06 


100 


11 


3501 


- 


51 




II 


3501 










| 


50 


STO 5 


3305 


10C 


II 


3501 



NOTE: Before putting a program into memory, after switching to W/PGM mode, press f 

and then PRGM to clear memory. 



CARD 3 



CARD 4 



Compute time-of- flight 



Key 

Entry 

STO 5 
X 

1 

CHS 

f 



Key 


Code 


Entry 


Shown 


RCL 4 


3404 


+ 


61 


STO 3 


3303 


R/S 


84 


LBL 


23 


B 


12 


STO 2 


3302 


R/S 


84 


LBL 


23 


C 


13 


STO 7 


3307 


R/S 


84 


LBL 


23 


D 


14 


g 


35 


RAD 


42 


RCL 5 


3405 


X 


71 


1 


01 


- 


51 


CHS 


42 


RCL 6 


3406 


-1 


81 


f ± 


32 


COS 


05 


RTN 


24 


LBL 


23 


E 


15 


g 


35 


DEG 


41 


RCL 7 


3407 


f 


31 


COS 


05 


RCL 2 


3402 


X 


71 


RCL 3 


3403 


X 


71 


ENTER 


41 


X 


71 


RCL 1 


3401 




81 


STO 8 


3308 


2 


02 


RCL 3 


3403 


T 


81 


RCL 2 


3402 


ENTER 


41 


X 


71 


RCL 1 


3401 




81 



60 



70 



80 



90 



/ 

STO 6 
g 

7T 

ENTER 

+ 

RCL 4 
D 

ENTER 

f 

SIN 
RCL 6 
X 



RCL 3 
D 

ENTER 

f 

SIN 
RCL 6 
X 



RCL 5 

RCL 1 
RCL 5 
X 
f 



R/S 

GTO 

E 

g NOP 
g NOP 
g NOP 



100 



Code 

Shown 

51 

3305 
71 
01 
51 
42 
31 
09 

3306 
35 
02 
41 
61 

3404 
14 
51 
41 
31 
04 

3406 

71 

51 

3403 

14 
41 
31 
04 

3406 

71 

51 

51 

3405 
81 

3401 

3405 

71 

31 

09 

81 

84 

22 

15 
3501 
3501 
3501 
3501 
3501 
3501 
3501 
3501 



Compute maximizing flight path angle 6 



10 



20 



30 



40 



50 



Key 

Entry 

2 

RCL 2 
ENTER 
X 

STO 5 
RCL 3 
X 

RCL 1 



STO 6 

RCL 4 
g 

1/x 

RCL 3 
g 

1/x 

RCL 5 

RCL 1 
X 

ENTER 

+ 

1 

+ 

STO 5 
g 

DEG 

0 

STO 7 
5 

STO 8 

LBL 

0 

RCL 5 

RCL 3 

RCL 4 

RCL 6 
X 

RCL 7 
f 

COS 

X 

gLSTx 

X 

RCL 7 
f 



Code 

Shown 

02 

3402 
41 
71 

3305 

3403 
71 

3401 

81 

51 

3306 

3404 
35 
04 

3403 

35 

04 
51 

3405 
81 

3401 

71 

41 

61 

01 

61 

3305 

35 

41 

00 

3307 

05 

3308 
23 
00 

3405 

3403 

3404 
81 

3406 
71 

3407 
31 
05 
71 

3500 

71 

51 

3407 

31 



60 



70 



80 



90 



100 



Key 

Entry 

SIN 

X 

RCL 5 
RCL 3 
RCL 4 

RCL 7 
f 

COS 

X 

ENTER 

X 



rr 

i 

+ 

RCL 6 



COS 

f 

TF 2 
R/S 
g NOP 
RCL 7 
gx<— >y 
STO 7 
0 

g x=y 
GTO 
0 
+ 

+ 

2 

STO 7 
g 

DSZ 

GTO 

0 

R/S 

5 

STO 8 
GTO 
0 

g NOP 



opt 

Code 

Shown 

04 
71 

3405 

3403 

3404 
81 

3407 

31 

05 
71 
41 
71 
51 
31 
09 
81 
01 
61 

3406 
81 

31 
09 

32 

05 

31 

81 

84 

3501 

3407 
3507 
3307 

00 

3523 

22 

00 

61 

61 

02 

81 

3307 
35 

83 
22 
00 

84 
05 

3308 
22 
00 

3501 



NOTE: Before putting a program into memory, after switching to W/PGM mode, press f 

and then PRGM to clear memory. 
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VII. The Appendix . 



As indicated in section HID, the maximum range for a 

specified cutoff altitude h„ and cutoff velocity v T is 

b L 

found by the following unconstrained optimization: 



S max “ S(6 opt> = max S(6) 



0<gs 



IT 



Using the equations in section IIIA, we want to find 

8 . such that 

opt 



3S 

86 



r E 8 8 (f T 



6 = 6 



= 0 



opt 



6 = 3 



opt 



We find that 



86 (f T f L } 



where 






II 

1 


lR , 
86 


1 t 1 


/p_ _ A 


er . o 3 

l 2ae 


K ) 
_ 


E T~ f L With 


respect 



r i=T,L 



vanish at 6=6 we must have either: 

opt 



1. 



^ y t y l } 6=6 . =0 

opt 



15 



or 



2 . 



3P 

33 



o 



3 = 3 



opt 



We find that 




2 2 . 
r v sin 

■Li Li 

y 



23 



0 



7T 

for 3 = 0, ± /2, ± ir, etc. Since we require that the 

launch point and the target point be separated by the line 

7T 

of apses, we require that 6 < 8^2 • This requirement is 
imposed because we must properly interface the disconti- 
nuity caused by the principal angle in the arc cosine used 
in steps 5 and 6 of section III A. Thus 3 = U /2 is the 
only angle of interest for which 3P/33 = 0, and this 
angle clearly gives rise to a minimum in range since the 
trajectory of the missile is vertical (straight up, and 
then straight down) . 

The value of 3 which maximizes S is found from the 
requirement that 



( Ym"Y T ) p _ o =0 
T L 3 ~ 3 opt 



This equation has defied attempts to find an analytic 
solution, but it can be rearranged and displayed as 
equation 1 in section III D, where the solution procedure 
is detailed. 
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